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INTRODUCTION 


The first scientific working group meeting of the MSFC Airborne 
Doppler Lidar Wind Velocity Measuranent Program was held August 25-26, 
1980, at the Marshall Space Flight Center. The working group meeting 
followed naturally from an exploratory meeting on the MSFC Airborne 
Doppler Lidar system held at MSFC in April 1980. This effort is part 
of the objective of the Severe Storms and Local Weather Research Program 
of the NASA Office of Space and Terrestrial Applications (OSTA). 

The purpose of the first scientific working group meeting was 
fourfold: (1) to identify flight test options for engineering verifica- 

tion of the MSFC Doppler Lidar; (2) to identify flight test options for 
gathering data for scientific/technology applications; (3) to identify 
additional support equipment needed on the CV-990 aircraft for the flight 
tests; and (4) to identify postflight data processing and data sets re- 
quirenents. The working group identified approximately ten flight 
options for gathering data on atmospheric dynamics processes, including 
turbulence, valley breezes, and thunderstorm cloud anvil and cold air 
outflow dynamics. These test options will be used as a basis for planning 
the fiscal year 1981 tests of the Doppler Lidar system. 

In addition to the minutes of the first scientific working group 
meeting, this report contains summaries of the material presented at the 
meeting, reports by the working group technical teams, and a meeting 
attendance list. 



MINUTES OF THE FIRST SCIENTIFIC WORKING GROUP MEETING OF 
THE AIRBORNE DOPPLER LIDAR WIND VELOCITY MEASUREMENT PROGRAM 


The first scientific working group meeting of the MSFC Airborne 
Doppler Lidar Wind Velocity Measurement Program was held on August 25-26, 
1980, at the NASA Marshall Space Flight Center. 

The morning of the first day was devoted to providing the working 
group members with information about the lidar system, its status and 
schedule, data processing, wind retrieval from Doppler lidar observations, 
an overview of the Cooperative Convection Precipitation Experiment (CCOPE), 
and CV-990 flight operations. During the afternoon of the first day and 
the early morning of the second day the working group divided into teams 
to: (1) identify and develop flight test options, (2) identify postflight 
data processing and data set requirements, (3) review real-time data display 
requirements, and (4) identify additional support equipment needed on the 
CV-990 aircraft for the FY81 flight tests. The results of these team 
activities were reported in the plenary session of the working group 
during the late morning of the second day. 

The convective phenomena team was chaired by Dr. James Telford of 
the University of Nevada; the secondary flows, boundary layer, turbulence, 
and wave team was chaired by Dr. Daniel Fitzjarrald of the Florida State 
University Geophysical Fluid Dynamics Laboratory. A list of the team 
members is provided as part of Appendix A. A float team consisting of 
representatives of the Doppler lidar hardware, data processing, and flight 
operations aspects of the program supported the activities of the science 
teams. 


The agenda of the meeting is included with the vugraph material 
presented by Dr. George H. Fichtl (see Appendix A). This material also 
includes a summary of the key results of the April 1, 1980, exploratory 
meeting on the scientific and technological applications of the MSFC 
Doppler Lidar System, the scientific opportunities to apply the lidar 
system, an updated science working group schedule, experiments, requirements, 
and a list of the working teams and their activities. Also included is a 
chart on the Doppler lidar characteristics which depicts how wind velocity 
data are acquired by the Doppler lidar. 

Mr. James Bilbro gave a status review of the Doppler lidar system 
and schedule. The basic points of the presentation were: (1) the status 
of the Doppler lidar hardware and modifications, (2) placement of the 
lidar and associated hardware on board the CV-990 aircraft, (3) various 
procedures for pulsing the laser beam, (4) data reduction options and 
error analyses in regard to pointing the laser beam, (5) range capability 
(i.e., out to approximately 10 km), etc. 

Dr. Harold Jeffreys and Mr. Carl Buck provided information on the 
real-time data processing and display of data. Also, the postflight data 
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processing issue was discussed (see Appendix B). The data display topics 
included the on-board aircraft data acquisition system (ADDAS) as well as 
the inertial navigation system (INS) and ancillary data (e,g., temperature/ 
dew point) to be obtained. The working group members expressed concern that 
the current concept for displaying the maximum intensity of the Doppler 
lidar backscatter signal might be difficult to interpret during flight 
operations. It was suggested that a separate CRT display should be used to 
display such information in view of the fact that new data are important 
relative to assessing Doppler lidar performance and for making scientific 
and operational decisions in real time. James Bilbro pointed out that a 
certain amount of modification of the currently planned on-board processing 
and data display is possible, so that a separate display of Doppler lidar 
backscatter intensity may be possible. The only constraints for imple- 
menting the separate display of Doppler lidar backscatter intensity are 
computer memory size and cost. 

Questions arose about the amount of time required to hardcopy the 
Doppler lidar data aboard the CV-990 aircraft after it is displayed on 
the output CRT (which also serves as the storage device of processed 
Doppler lidar data for the on-board data processing). It was pointed out 
that it takes 10 seconds to obtain a hardcopy of data, so that very little 
data would be lost during the copy routine. 

The postflight processing of wind components into wind vectors was 
discussed. It was pointed out that the forward and rearward looking 
Doppler lidar beams are separated in time, with the time differential 
increasing with range. It was further noted that combination of the 
forward and rearward looking beams into wind vectors is a matter of 
scientific interpretation. Therefore, the working group suggested that 
data sets of forward and rearward looking Doppler lidar winds should be 
provided, in addition to a data set in which the forward and rearward 
looking beams are combined in a rational manner to yield total wind vectors. 

In the context of defining cloud boundaries, questions arose 
(Dr. Richard Doviak and Dr. William Vaughan) concerning the distance the 
laser beam could penetrate into a cloud. James Bilbro said that in most 
convective cloud cases the Doppler lidar beam will have negligible pene- 
tration; however, if clouds are significantly thin (e.g., cirrus), signifi- 
cant penetration can occur. The answer will depend on the drop or crystal 
size distribution and liquid water or ice content within the cloud. It 
was also noted by James Bilbro that, for a typical convective cloud, a very 
pronounced change in backscatter signal intensity will occur and that this 
return could be used to define cloud backscatter and cloud boundary in 
the postflight data processing activity. 

Dr. Robert W. Lee provided a presentation on wind retrieval from 
Doppler lidar observations. His vugraphs are provided in Appendix C. 

This presentation clearly pointed out the need for precise specification 
of postflight data set requirements relative to accuracies, confidence 
bands, combination of forward and rearward looking wind components to 
obtain total vector winds, and calculation of kinematic quantities, such 
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as vorticity, divergence, etc. Appendix D is a report by Robert Lee on 
wind field retrieval from Doppler lidar observations which includes much 
of the material presented by Dr. Lee at the meeting. 

Dr. Lavon J. Miller provided an overview of the CCOPE. Appendix E 
is the material used in his presentation. Extensive discussion took place 
after this presentation relative to (1) flight operations and (2) the kinds 
of data that could be acquired with the Doppler lidar at the CCOPE. It 
is clear from the discussions that close coordination and extensive planning 
with the CCOPE on the part of the Doppler lidar team members will be 
required if the Doppler lidar is used at the CCOPE. In this regard Lavon 
Miller requested that representatives from the Doppler lidar team partici- 
pate in the October CCOPE aircraft operations planning meeting to be held 
at NCAR, Boulder, Colorado. Current plans include the attendance of 
George Fichtl and John Kaufman at this meeting. Furthermore, during the 
CCOPE, close coordination of the CV-990 flight operations with the CCOPE 
operations center in Miles City, Montana, will be required. This could 
prove to be problematic because the CV-990 will most likely be stationed 
at Rapid City, South Dakota. However, it is not a problem that cannot be 
solved. 

The CCOPE offers the opportunity to acquire ground truth wind velocity 
measurements with the Doppler radars in the CCOPE radar network for verifi- 
cation of the Doppler lidar system. In addition, scientific opportunities 
relative to performing fundamental studies with the Doppler lidar may be 
possible at the CCOPE in regard to gust front dynamics, anvil dynamics and 
precipitation efficiency, dry air entrainment associated with cumulus 
turrets, and feeder flow under small cumulus clouds, as well as boundary 
layer studies in prestorm conditions. This is discussed in more detail in 
the report by the convective phenomena team (Appendix G). 

George Alger commented on how CV-990 operations will take place. He 
pointed out that the CV-990 has approximately a 7-hour flight endurance. 
Typically for a flight test a 1-hour period will be required for flight 
planning, and the flight would take place approximately 2 hours after the 
planning phase. Thus, weather forecasts will play a key role in our flight 
activities. George Alger pointed out that, once the CV-990 is in the air, 
changes in the flight plan are possible. However, such changes in flight 
plans must be coordinated with James Bilbro and George Alger to coordinate 
the changes with the pilots. Mr. Alger stated that the first two to three 
flights will most likely be used for getting "bugs" out of the instruments 
and familiarizing the flight crew with the operation of the Doppler lidar. 

After these discussions, George Fichtl gave a presentation in prepara- 
tion for the team meetings that took place in the afternoon and the next 
morning (August 26). His presentation charts are given in Appendix F. 

The reports of the team meetings are provided in Appendices H through 0. 
Reports by Dr. Richard Doviak (NSSL) and Dr. Randall Koenig (NSF) provide 
details on recommended flight test options, recommended flight paths, 
on-board data displays, postflight data processing, scientific and technical 
rationale, and purposes (Appendices I and J. 
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The meeting proved to be successful in accomplishing the goals set 
forth earlier in these minutes. The identified flight options and other 
inputs will be used to prepare a flight option document. This document 
will summarize the experiments and will be used: (1) for determination of 
Doppler lidar operational settings, e.g., pulse length, number of pulses 
to be used in calculating averages, etc., and (2) for CV-990 flight 
planning. 


Questions regerdinq 
(205 453-0875) or John W. 




these minutes may be directed to George H. 
Kaufman (205 453-3104). 


Fichtl 


I 
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APPENDIX A 


PART I OF PRESENTATION BY 
DR. GEORGE H. FICHTL, MSFC 
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DR. ROBERT W. LEE, LASSEN RESEARCH 
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Abstract 


An approach to the retrieval of a vector wind field from 
Doppler lidar observations is developed in general terms. The 
field of radial velocity measurements from each look angle is 
modeled by a smooth surface, the parameters of the model being 
determined from the data by least-squares techniques. The 
vector wind field and higher-order fields are obtained from the 
two modeled surfaces. Estimated measurement errors are taken 
into account, and error estimates are available for all output 
data sets. 


D-2 



Table of contents 


I, Introduction 1 

II* Problem definition 2 

A, Problems with the data 2 

B. Desired results 2 

C* The solution in general terms 2 

D. Steps required for implementation 3 

III. Expanded problem definition 4 

A* Specific nature of data 4 

B. Inversion goals 6 

C* Generalized solution 7 

D, Solution implementation 8 

IV. Detailed Implementation 9 

A* Data editing 9 

B* Data smoothing 10 

1* General 10 

2. Continuous surfaces 10 

3. Surface adaptability 11 

4. Suitable models 11 

5. Model fitting 12 

V. Short topics 13 

A. Spectral width an<l amplitude 13 

B. An alternate solution 13 

VI. Algorithm Implementation 14 

A. Implementation 14 

B, Evaluation 14 

VI I . Conclusion 16 


D-3 



List of figures 


1. Probability distribution of Doppler estimator error 5 

2. Exaggerated observation geometry 5 


D-4 


I • Introduction 




This report is concerned with the problem of retrieving vector wind 
measurements in a plane from radial wind measurements made in that plane 
using two different look angles. There are several possible approaches to 
this problem, but certain characteristics of the data rCQulre a certain amount 
of cafe in the selection of the approach. 

The measurements are made by an airborne Doppler lidar system. The 
sensitivity of this system is limited, and the windfleld tracer - naturally 
occurring aerosols - is not always present in sufficient quantities for a 
satisfactory return. For this reason the radial velocity measurements will 
vary greatly in accuracy, causing errors in the inferred vector field which 
will be magnified by the less than perfect orthogonality of the two look 
angles. Finally, the utility of the higher-order attributes to be derived 
from the vector windfleld - vorticity and divergence, for example - will be 
limited by noise and error in the Inversion process. 

These considerations clearly suggest an inversion algorithm which is 
tolerant of gross measurement errors and which minimizes the effects of random 
errors in the data. Achieving these goals without greatly reducing the 
resolution of the measurements requires care. 

Section II below consists of a more precise definition of the inversion 
problem. Elements of this definition are presented in detail in section III. 
The primary inversion steps of editing and smoothing are considered in section 
IV. Two related topics are discussed briefly in section V, while the steps 
required for implementation and evaluation of the suggested algorithm are 
outlined in section VI. This final step - the evaluation of the solution 
algorithm - is just as Important as the selection of the algorithm itself. 

Only when the error characteristics of the solution are known can one 
interpret the measurements with confidence. 
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II* Problem definition 


Simply stated, the goal is to derive a two-dimensional vector field from 
two scalar fields. In practice, most of the data manipulation will take place 
on the separate scalar fields (see section V.B, however, for an alternate 
strategy), and conversion to a vector field becomes a final, trivial operation. 

Since a solution to the problem posed depends upon the definition of that 
problem, some care must be exercised in characterizing the problem itself. Of 
particular importance here are certain problems inherent in the measurements, 
and the need for a considered definition of the properties desired in the 
solution. 

II. A. Problems with the data 

The measurements of the two radial velocity fields may be sparse (missing 
data points), and they certainly will vary greatly in quality. These two 
attributes will be pre-eminent in the selection of a solution algorithm. In 
addition, the measurements are made on irregular grids, with no agreement 
between the measurement points in the two scalar fields. Finally, the two 
scalar observations are not orthogonal, making the conditioning of the scalar- 
vector transformation less than ideal. 

II. B. Desired results 

It is not easy to quantify the characteristics required in the solution 
data set, but it is important to attempt to do so. The solution algorithm 
will be optimized according to these criteria, so obviously the criteria must 
be appropriate. 

The general requirement is for smooth flow fields of known accuracy. 

These fields must be useful in a visual sense, and the errors in the statistical 
properties of these fields must be known and acceptably small. 

II. C. The solution in general terms 

Clearly the path to the solution is one of data smoothing and Interpolation. 
These steps must be accomplished in a way that is optimal given the data 
characteristics and the desired goals of the solution field. In addition, the 
solution must be efficient in that it makes use of all a priori information 
about the measurements. 
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II. D. Steps required for solution implementation 

In addition to careful statement of the problem and the desired solution 
goals, selection of an algorithm requires a procedure for evaluation of the 
algorithm according to those goals. Such an evaluation will require operation 
on real and simulated data sets, with visual and statistical interpretation of 
the end products. Under these conditions the algorithm can be ”tuned" and 
optimized for the actual characteristics expected of the raw data. 
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Ill* Expanded problem definition 

III.A. Specific nature of data relevant to optimal field retrieval 

As mentioned above> the data is sparse* This may be due to range 
attenuation, lack of aerosols, or both. While missing data points will be a 
particular problem at long range, they can occur at any range. 

The quality of the measurements will vary greatly* The Doppler estimator 
used for measurement of radial velocity has certain (known) error characteristics 
depending upon the slgnal-to-noise ratio of the signal return and the spectral 
width of the signal. The probability distribution of this error is sketched 
in figure 1. It consists of two parts: a normally-distributed terra N and a 
uniformly-distributed term U* At very low signal-to-noise ratios the uniformly- 
distributed term gives rise to *’wild** estimates of radial velocity, which 
account for most of the error variance. 

In addition to estimation errors produced by the Doppler processor there 
is an Inherent sampling error present regardless of the signal-to-nolse ratio. 
Since the Doppler return has finite bandwidth any finite realization of that 
return will be subject to an error variance due to sampling. At high signal- 
to-noise ratios this is the dominant error term. 

The locations at which radial velocities are observed do not form a 
regular grid. The angles at which the observations are taken will vary 
slightly due to aircraft dynamics, as shown in exaggerated form in figure 2. 

The gridpolnts of one scalar field will bear little or no relationship to 
those of the other field* The observation angles will not differ by the desired 
90 degrees, but by approximately 40 degrees* This lack of orthogonality will 
magnify data errors in the vector component parallel to the aircraft track* 

The spatial sampling of the laser beam presents further problems* The 
beam resolution cell is long (~300 m) and narrow (--30 cm radius); each 
measurement will consist of the average of a number of such cells displaced 
horizontally by approximately 1 meter. Such a resolution volume will average 
effectively in one spatial dimension (range), but not in the other two* The 
resulting under-sampling of the spatial windfield (particularly in the vertical 
dimension) will lead to an increased error variance in the data from aliased 
energy in the spatial windfield* This is further complicated by the fact that 
the measurements are not necessarily made in the desired plane, but at a 
horizontal angle subject to random excursions about zero* Indeed, the effect 
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P(v) 



Flgvre 1: Probability distribution of Doppler estimation 

errors for a certain slgnal-to-noise ratio. N and U are 
normally- and uniformly-distributed variance components. 



Figure 2: Exaggerated observation grid. X and 0 represent 

backward and forward look angle data points. 
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will be reduced by intentionally jittering the vertical axis of the system to 
achieve some averaging in that dimension. 

If the turbulent parameters of the wlndfield can be estimated - and this 
should be possible from the data itself - then the variance due to under- 
sampling can be estimated. 

Finally, there is a d-ata problem due to the fact that the two scalar 
fields do not result from simultaneous observations, /rhe time lag between the 
observations varies with range, being on the order of 30 seconds at maximum 
range. This is sufficient time for the wlndfield to translate by one or more 
resolution cells, and it may be necessary to advect the measurements according 
to the local mean v/ind vector for an appropriate time interval (which may vary 
with range) in order to achieve the required time registration of the 
measurements. 

I I I . B. Inversion goa Is 

The first use of the vector flow field will be visual, in the understanding 
of flow fields in the clear-air vicinity of severe storms. Thus the inferred 
vector field must be sufficiently smooth to allow the eye to trace parcel 
flow. Stated another way, the continuity between neighboring vector estimates 
must be reasonably high. 

Aside from visual aesthetics the vector field must possess certain 
statistical properties. Useful vorticlty and divergence fields must bo 
obtainable from it. This in turn requires that the errors in the estimated 
divergence and vorticlty be below some threshold (perhaps a certain percentage 
error) for a substantial percentage of the observations, and in addition that 
the errors themselves can be estimated. Obviously a vorticlty estimate of 
O.Ol/scc is not useful if the standard error is O.l/sec, or if the error is 
unknown. The allowable errors In vector wind and higher-order functions, 
expressed as percentages at some confidence level, are very important inputs 
to the algorithm optimization and evaluation process. 

Finally, the inverted data must be as complete as possible. Gaps in 
the derived fields must be filled where possible by Interpolation, even though 
this may reduce resolution in the vicinity of the missing data points. Caro 
must be exercised in such interpolation to ensure that good data points aro 
not contaminated by bad. 
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III.C, Generalized solution 


In general smoothing and interpolation are achieved through the reduction 
In the degrees of freedom of the solution field below those of the input data, 
and through careful selection of those degrees of freedom retained in the 
solution. Only by sacrificing degrees of freedom - and here resolution - In 
the solution field can data reliability be enhanced and error variance estimated. 
Stated in spatial terms, resolution is reduced to obtain greater data stability. 
In spectral terms, knowledge of field components of high spatial frequency is 
sacrificed so that knowledge of lower spatial-frequency components will be 
enhanced. Clearly a balance must be struck somewhere in between the extremes: 
perfect knowledge of the mean wind on a 10x10 km square is of little value in 
the severe-storm program, as is no knowledge on a 100x100 m grid. 

The general process of smoothing and interpolation is one of modeling. 

One selects a mathematical framework by which to model the output field, and 
determines the parameters of that model from the measurements. 

The success of the smoothing process depends to a large extent upon the 
suitability of the model selected. The model must be appropriate in several 
respects. First it must be able to represent the natural features of the 
windfield adequately. It must be possible to control the spatial-frequency 
response of the model readily, to allow control over the smoothness of the 
solution. Finally, it must be mathematically tractable: a model is not useful 
if it takes an hour of computer time to invert a minute’s worth of data. 

Once a model is selected the parameters of that model may be determined 
from the data. Such a parameter fit may be achieved in a straightforward 
manner using least-squares techniques, but note that this requires use of 
objective error criteria. The solution is optimal in terms of these criteria, 
but one must make sure that the criteria are appropriate. In general these 
criteria will involve some compromise between spatial resolution and data 
stability, between smoothness or continuity and error variance. 

Given that a model has been selected and the solution obtained for a 
certain set of error criteria, it remains to evaluate the probable errors in 
the various end products - wind vectors, vorticity, etc. 
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III.D. Solution implementation 

The generalized solution outlined above can be implemented in the following 
sequence: 

1) Obtain data. In addition to recording the raw data, this step 
Includes obtaining all ancillary data which will be useful in data interpretation 
time and location, look angles, other meteorological data, etc. 

2) Establish data reliability. From the slgnal-to-nolse ratio estimate 
and other parameters, estimate the probable error of the velocity estimate. 

This error may consist of two parts with different probability distributions. 

3) Editing. The data must be put into a form that the smoothing algorithm 
can use. In addition to assigning a weight to the data point reflecting its 
probable error, and assigning to that point coordinates, spectral width, etc., 
an additional operation is desirable. The measurement may be compared with 
neighboring points (in two dimensions) as a test of measurement continuity. 

If the measurement is discrepant its weight may be reduced. This process will 
remove to a great extent the effects of "wild" estimates produced by the 
uniform portion of the Doppler estimator error. This operation is explored 
in more detail in section IV. A. 

4) Smoothing. This step includes solving for the model parameters in 
terms of the weighted data, interpolating where required, and reducing the 
resolution of the measurements where data quality is low. Estimates of error 
variances should be carried through this process. Finally, the solution field 
can be sampled on any desired grid. This operation is described in detail in 
section IV. B. 

3) Produce vector field. The two scalar fields may be combined to form 
the vector flow field, again carrying through the estimated errors. 

6) Produce end products. Hlgher~order fields may be obtained through 
operations upon the vector flow fields, in each case carrying error estimates 
through the process. 

7) Evaluation. The last step is to evaluate the utility of the resulting 
end products. If there are serious problems with them, it must be determined 
from error propagation which aspect of the raw data most seriously compromises 
the result. If that data aspect cannot be corrected, it should be determined 
whether or not alterations in the model can reduce the effects. 


D»12 


IV. Detailed Implementation 
IV. A, Data editing 

The goal of the data editing process is to produce data of known error 
characteristics for the smoothing algorithm. All Information available must 
be brought to bear in evaluating a given data point. The following list 
Includes the most important factors. 

1) Slgnal-to-noise ratio. From the signal amplitude estimate at the 
range gate of interest, in comparison with the amplitude estimates at very 
large ranges (where no signal is expected), an estimate of the slgnal-to-nolse 
ratio can be obtained. This estimate can be used in conjunction with the 
(known) error , characteristics of the Doppler estimator to produce a probable 
velocity error estimate consisting of two parts as suggested above: a normally- 
distributed component and a uniformly-distributed component, 

2) Spectral width. The Doppler estimator produces as a matter of course 
an estimate of the signal spectral width. Spectral width enters into the 
Doppler processor error equations. Note however that useful estimates of 
spectral width are not produced at very low signal-to-nolse ratios. 

3) Continuity, Continuity may be used in two dimensions as a check 
upon data consistency. Continuity tests may be applied to amplitude and width 
estimates as well. In a typical case the weighted median value of the eight 
neighboring points might be compared with the point in question. Note that 
the median or most probable value is more useful here than the mean value, 
since the mean can be severely disturbed by a single bad data point. 

4) Constraints. The expected characteristics of the wlndfield can be 
used as a further check upon data integrity. For example, a constraint upon 
velocity gradients (shear) can be used as an input for continuity tests. 

Limits may be set upon maximum values of velocity as a test for reasonableness. 
As with all constraints, care must be exercised to ensure that actual features 
of the wind field which were not expected are not obscured. Use of adaptive 
or interactive constraints can achieve this goal. 

Should a given data point fail one or more tests for reliability, the 
weight of that point may be reduced, or in severe cases a missing data point 
may be declared. 

Note that the editing process can be combined with the smoothing or 
filtering process. A first fit of the data points gives a trial solution and 
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a deviation for each data point. These deviations are a measure of continuity 
and can be used to alter the weighting given the data. A second Iteration of 
the solution gives a revised output field. 

IV. B. Data smoothing 
IV.B.l General 

Data points can be considered in isolation, but since smoothing Implies 
that each data point has an Influence upon its neighborhood it is useful to 
consider the two-dimensional measurements as forming the height of a two- 
dimensional surface. The process of data smoothing then becomes one of 
fitting a surface of a certain character as nearly as possible to the 
measurements, in* (for example) a least-squares sense. 

Each form of data manipulation can be interpreted in terms of a surface 
of a certain type. For example, point data may be considered to form a 
surface composed of blocks centered at the measurement points, the height of 
each block indicating the value of the measurement at that point. That is, 
the data point is the altitude of the surface for that grid square. A 
continuous surface can be created by placing the data points at the vertices 
of the surface, with straight lines Joining the vertices defining the surface 
(that is, linear interpolation between the data points, with grid rectangles 
broken into two triangles by a single diagonal) . Surfaces formed with 
continuous first- or higher-order derivatives require the overlapping influence 
of several measurement points at each point on the surface. 

IV. B. 2 Continuous surfaces 

Continuous surfaces may be modeled by many analytic or elementary 
functions. The most commonly used functions are polynomials (including splines, 
Hermlte and other orthogonal polynomials), Fourier series, Bessel functions 
and spheroidal functions. The choice of a basis function depends upon: 

1) Suitability for the problem. Some functions lend themselves to a 
particular coordinate system or situation. For example, Bessel functions are 
often appropriate for cylindrical coordinates, and Fourier series for band- 
limited functions. 

2) Mathematical ease of use. Polynomials, for example, offer few 
difficulties in integration or differentiation, no convergence problems, etc. 
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3) Parameter flexibility. The degrees of freedom of some functions 
can be easily ’’tuned*’ to control the parameters of the surface. For example, 
the small-scale wiggles of a surface modeled by a Fourier series are easily 
controlled by limiting the order of the series. 

IV. B. 3. Surface adaptability 

Since the quality of the data varies from point to point on the surface, 
it may be reasonable to allow the nature of the surface to vary as well. That 
is, in regions of high data quality the smoothness constraint on the surface 
may be relaxed to allow the model to represent smaller-scale features. 

Conversely, in regions of poor or missing data surface smoothness must be 
constrained e^en further to preserve surface continuity. This trade-off 
between surface smoothness and resolution may be made in an adaptive manner, 
with the algorithm Itself sensing the need for constrained smoothness. 

IV, B, 4, Suitable models for flow fields 

Due to its easily-controlled spatial-frequency response, a Fourier 
surface is attractive. However, the difficulty of incorporating data of 
varying quality, on a non-uniform grid, is substantial* The model of choice 
is a locally-defined polynomial with a basis function of limited extent. 

Such a model offers ease of solution using least-squares techniques, no grid 
problems, controllable basis size and smoothness, and continuity to any 
desired derivative. Suitable basis functions would include linear, quadratic 
or higher functions over limited (sliding) basis regions, or spline functions* 

As an example, consider the lowest-order surface fit, A region of 
influence is defined around a point for which surface height is to be estimated. 
Data points in this region of influence are summed in a weighted average, the 
weights being derived from the error variances assigned to those points, with 
(for example) an additional weighting function formed by a two-dimensional 
Gaussian centered at the estimation point. This weighted averaging is equivalent 
to fitting a local plane to the data in the vicinity of the estimation point. 

To achieve the accuracy desired, the size of the region of Influence (defined 
by the two-dimensional Gaussian weighting function) can expand or contract as 
required to enclose a suitable number of measurement points. Such an approach 
is easily Implemented, and sliding the Gaussian region of influence around the 
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plane gives a continuous estimation surface. However, with this simple approach 
shear in the windfield cannot be fitted by the model at each point; the result 
is a poor fit requiring increased smoothing. 

Use of the next higher order surface solves the shear problem: at each 
estimation point one fits the height and slope of a plane surface. Shear is 
no longer « problem - the fitting errors are limited to curvature and higher- 
order derivatives. 

At some point increasing the order of the model, (that is, increasing the 
degrees of freedom in the solution) Increases noise in the surface beyond a 
tolerable level. The optimum surface order remains to be determined; there 
will be a compromise between higher order auid reduced region of influence which 
must be determined by simulation. 

Splines are a particularly attractive form of polynomial basis function 
since the approximating functions are easily constrained to be continuous on 
the boundaries between grid points. 

IV. B. 5 Model fitting 

Once a model has been selected and the controllable parameters defined, 
it remains to determine those parameters. The most suitable solution technique 
is the linear least-squares approach. The variable weights of the data points 
are easily taken into account, along with additionally-imposed geometrical 
weighting. One particular advantage of this approach is the availability 
with the solution of an estimate of the solution error variance. 

Once the parameters of the surface have been determined, that surface 
may be sampled at any desired grid pattern. 


V* Short topics 

V,A. Spectral width and signal amplitude 

Although they have received little attention thus far, the signal amplitude 
and spectral width are also measured by the Doppler processor. These two 
quantities may also be considered to form solution surfaces, and the same 
techniques described above may be applied to the estimation of the parameters 
of these surfaces: editing, smoothing and interpolation. 

Additional redundancy is present in these measurements, since information 
from the two look angles may be combined. 

Note that a portion of the apparent spectral width may be contributed by 
horizontal velocity shear within the target volume. Since the velocity field 
is being determined independently, it is possible to correct for this con- 
tribution under the assumption that the shear variance and the spectral width 
add incoherently. 

V.B. An alternate solution strategy 

While this report has treated the data from the two look angles as being 
independent until their combination in the vector field, another approach is 
possible. 

The wind field model may assume a single surface as a potential field. 

The measurements become directional derivative estimates of this surface, and 
techniques for surface reconstruction from derivative information can be used. 

Note however that this process is strictly valid only when the divergence 
of the actual wind field is zero. Thus the potential field so derived will 
naturally produce a zero-divergence field. Divergence may then be recovered 
from the measurements by a second-stage solution, solving for a divergence field 
from the difference between the measurements and the inferred zero-divergence 
field. There may be a problem here since the divergence and circulation may be 
locally correlated. 



VI. Algorithm implementation and evaluation 
VI. A. Implementation of the surface-fitting model 

This section Is an outline of the steps required to take this solution 
technique from the generalized concept described In this report to the 
functional stage. The primary questions to be resolved at this point are: 

1) Definition of the most appropriate surface model 

2) Definition of the appropriate solution technique for that model 

These two questions are Inextricably Joined. Their solution will arise 

through an Iterative process of evaluation and optimization. 

Once a model and a solution technique are chosen, they must be "tuned up" 
with reference tp the practical problems of the data at hand. This must be 
done with the data end use firmly In mind, employing a well-defined evaluation 
technique and a set of evaluation criteria. Such tuning will determine the 
appropriate editing and weighting schemes, grid sizes and spatial resolution. 

The suitability of the resulting algorithms Is critically dependent upon 
the accuracy and realism of the evaluation technique, discussed In the next 
section. 

VI. B. Algorithm evaluation 

Since the solution technique will be optimized through Interaction with 
a process of evaluation, the technique will ultimately by optimal only In terms 
of that evaluation procedure. Only If the evaluation procedure reflects the 
realities of the data and the wind field can one expect the solution algorithm 
to be optimal for the data. 

In addition to providing a test bed for optimizing the solution algorithm, 
evaluation provides two Important byproducts: 

1) Confidence In the results. If the user can take a real or S 3 mthetlc 
wind field, probe It with a simulated Ildar system, contaminate the data with 
reasonable errors, and still retrieve a useful approximation of the original 
wind field, then he can have some confidence In using the algorithm upon data 
for which there Is no confirming data. 

2) Error propagation. By use of simulation the errors In user products 
can be estimated In terms of the errors In the raw data. User products without 
error variance estimates are of marginal utility; this Is especially true of 
higher-order products such as convergence. 
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The following items may be taken as defining the components of an 
evaluation program: 

1) Goals* A set of target goals should be established, In probabilistic 

terms* For example, one might desire that the vector wind components be 

-3 -1 

measured to 2 m/s 90% of the time, or that vortlclty be accurate to 10 s 
on a l->km scale* 

2) Input data sets. Both simulated and actual wind fields (taken from 
multiple-Doppler observations) are of value - the former for their controlled 
nature and the latter for their realism* Obviously the statistical properties 
of these fields must be accurately known, 

3) Signal-to-noise ratio* For simulation purposes, realistic slgnal-to- 
noise ratios must be assigned to the data points on a random basis* This would 
include range variation, dropouts, Rayleigh statistics, etc. This signal-to- 
noise ratio will be used to assign probabilistic errors to the radial velocity 
simulations, so it is important that they be realistic* 

4) Wild measurements* In the transition from signal-to-noise ratio to 
velocity error, an appropriate number of tot ally-random estimates must be 
included to reflect the component of Doppler estimator error which is 
uniformly distributed* 

5) Geometry* The grid points and look angles should be varied in a 
random way with reasonable values of variance. 

6) Solution* Given the velocity field as probed by the synthetic Ildar 
system - noise and all - an estimate of the original velocity field may be 
obtained by using the solution algorithm under test* 

7) User products* The output wind field estimate may be transformed 
into the desired end products: visual fields, statistics, higher order fields* 

8) Evaluation* The errors and utility of the user products must be 
assessed through comparison with the initial data set, using the evaluation 
goals as criteria for success. 

The results of this evaluation may suggest alterations in the model or 
solution technique, or may suggest that certain user products cannot be 
reliably obtained from data of the quality simulated. By varying the charac- 
teristics of the input data set, the sensitivity of the inversion process to 
data problems of a given type may be determined. These sensitivity factors 
may suggest certain constraints upon the experimental operation, in order to 
improve recovery of a given user product. 
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VII. Conclusion 


This report has suggested an approach to the retrieval of a wind field 
estimate from Ildar measurements. This approach seems likely to draw the 
maximum amount of useful Information from that data. Note however that some 
degradation of system resolution Is required. 

The emphasis of this approach Is upon error analysis at all stages of the 
solution. It Is felt that user products (wind fields, divergence fields, etc.) 
without explicit error estimates and confidence levels are of marginal value. 
This Is particularly true of smoothed fields. A smoothed random field cannot 
be distinguished from the smoothed fields reported by dual Doppler observers, 
and one should have no confidence In such highly mathematical products unless 
shown an error propagation example. 

With such error analysis techniques, Including a carefully planned 
evaluation technique, one can be confident that one will know when the derived 
wind field has significance. This apparently modest claim Is highly important 
when an experiment Is likely to have marginal results: it is far more preferable 
to have a few good wind fields of known reliability - even if they represent 
only a small portion of the measured fields - than to have results of doubtful 
validity for all the data sets. 
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APPENDIX E 


MATERIAL PRESENTED BY 

DR. LAVON J. MILLER, NCAR, BOULDER, COLORADO 



Cooperative Convective Precipitation Experiment 


•Colocated with Dept, of Interior/ Water and Power Resource 
SERVICE/ HIPLEX 


•SE Montana at Miles City 


•Improve understanding of the physics of convective 

PRECIPITATION 


•Major emphasis is placed on obtaining a good description 

OF WHOLE CONVECTIVE PRECIPITATION SYSTEM 


•Framework within which single/ significant/ tractable 

PROBLEMS ARE INVESTIGATED 
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CANAD A 




Scientific Objectives ! 


’Hydrometeor evolution 
•Precipitation efficiency 
•Origins of ice 
•Entrainment or mixing 
•Storm structure and environment 
•Storm initiation 
•Atmospheric chemistry 


•Storm electrification 




Hydrometeor Evolution; 


•Growth through ice process 


•Trajectories/ growth environments and size 
distributions of graupel/ rain and hail 


•Graupel as embryo of hail and rain 


•Ice process is too slow to form precipitation 
within moderate (UO ms"^) updrafts in 
straightforward WAY 


•Precipitation-sized ice particles are observed 
IN moderate to strong updrafts 


Precipitation efficiency: 


•Part of more general water budget specification 

•Amount of water vapor input versus amount of 
precipitation output 

•Major factors in moisture budget of clouds 

Upward flux of vapor in updraft 
Downward flux of precipitation in downdraft 
Evaporation within precipitation shaft 
Conversion of cloud droplets to precipitation 
Entrainment or mixing of clear air into cloudy air 
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Origins of ice: 


•Dominant mechanism of precipitation formation is 

DIFFUSIONAL GROWTH OF ICE CRYSTALS FOLLOWED BY 
ACCRETION (riming) OF CLOUD DROPLETS BY ICE 
PARTICLES 


•Environment within which early formation of ice 

OCCURS 


•Which ice particles become precipitation embryos 


•Relationship between in-cloud ice particle concen- 
trations AND ICE NUCLEUS CONCENTRATIONS FOUND 
IN THE INFLOWING AIR 
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Entrainment or mixing: 


•Drier environmental air mixes with cloudy air 
•Process occurs both laterally and vertically 
•Fully three-dimensional and time-varying 
•Scale sizes and partitioning of relative importance 
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Storm structure and the environment: 


•General description of physical characteristics 

OF STORMS AND THE ENVIRONMENT WITHIN WHICH THEY 
DEVELOP/ INTENSIFY AND DECAY 


•Encompass (realistically) scales from small convective 

■ ELEMENTS OR CELLS ('^ 2S KM) TO STORMWIDE ('^10-50 KM) 


•Dynamics of cell genesis - growth and interaction as 

WELL AS INTERNAL MECHANISMS 


•Processes acting to determine cellular structure 
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Storm initiation: 


•General problem of initiation of convection 


•Intensification and organization of broad scale 

CONVECTION INTO SPECIFIC STORMS 


•Roles of elements such as topography/ atmospheric 

WAVES/ EXISTING STORMS/ SYNOPTIC FRONTS AND 
TURBULENCE 
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Atmospheric chemistry: 


•Effect of clouds on chemistry of the atmosphere 


•Cloud AND precipitation chemistry in Great Plains 
compared to eastern United States 


Storm electrification: 


•Climatology relating lightning to stages of cloud 

DEVELOPMENT 


•Specification of electric field development above 

CLOUD TOP 


•Study major thunderstorm electrification mechanisms 
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Facilities: 


•Radar - Surveillance and 7 Doppler radars (2 are 
dual-wavelength) 

•Surface - WPRS '^100 mesometeorological stations/ 50 
AT ^0 and 20 KM SPACING/ NCAR (PAM) '^28 
stations plus remaining WPRS STATIONS AT 
-^8 KM spacing 


•Upper air - WPRS (2) and NCAR (2) sondes with NASA 
(10) AT '^50-60 KM spacing 


•Satellite - Colorado State University collects and 

ANALYZES VISIBLE AND INFRARED FROM GOES- 
EAST VISSR 


•Aircraft - 12 to 13 powered and 1 glider 
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Radars : 


•NCAR/Field Observing Facility 

CP-2 (S-BAND DoppleR/ X'-band incoherent) 
CP-3 (C-BAND Doppler) 

CP-A (C-BAND Doppler) 


•Univ. Chicago - Illinois State Water Survey 
CHILL (S-BAND DoppleR/ X-band incoherent) 


•NOAA/Wave Propagation Laboratory 

WPL-C (X-BAND Doppler) 

WPL-D (X-BAND Doppler) 

WPL-E (K-band DoppleR/ dual polarization) 


•HIPLEX 

SWR-75 (C-BAND incoherent) 
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DISTANCE EAST OF CP-2 (km) 



Aircraft: 


•National Center for Atmospheric Research (NCAR) 

Beechcraft Queen Air, N30AD 
Beechcraft Queen Air, N306D 
North American Rockwell Sabreliner, N307D 

•NCAR/National Oceanic and Atmospheric Administration (NOAA) 

ScHWEiZER 2-32 Sailplane, N9929J 
Cessna 180, N52032 

•Univeristy of Wyoming 

Beechcraft Super King Air 200, N2UW 
Beechcraft Queen Air, NIOUW 

•South Dakota School of Mines and Technology (SDSM&T) 

North American T-28, N510MH 

•National Research Council of Canada (NRC) 

DE Havilland Twin Otter, CF-POK-X 

•Colorado International Corporation (CIO 
Lear JET 23, N88TC 

•National Aeronautics and Space Administration (NASA/Goddard 
Space Flight Center 

WB-57F 

•University of North Dakota (UND) 

Cessna Citation II 

•Commercial Operator 

Photographic/Chaff aircraft 
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Aircraft measurements; 


•Position ~ multiple aircraft positioning system 
(' v-10 trackable) 

•PRESSURE/ TEMPERATURE AND DEW POINT 

•Winds - horizontal/ vertical and turbulent 

•Cloud and precipitation particles 
LIQUID water content 

•particle size distributions (particle measuring 

SYSTEM PROBES/ CAMERA/ FOIL IMPACTOR/ HAIL 
SPECTROMETER/ ICE PARTICLE COUNTER/ PARTICLE 
REPLICATOR) 

•Nuclei populations 

MEMBRANE FILTERS/ IMPACTOR SLIDES/ AND EXPANSION 

chamber 

•Cloud photographs 

•Surface infrared temperature 

•Electric field and electrical conductivity 

•Cloud top temperature and altitude 

•Integrated cloud liquid water content and precipitation 

PHASE 


Experimental approach: 


•Study entire lifetimes of convective clouds 

•Field operations organized into three (arbitrary) 
stages 

Prestorm 
Early storm 
Mature storm 
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SFC 


Prestorm: 


•Boundary layer heating and moisture convergence that' 

START CONVECTION 


•Mechanisms of continuance or regeneration of storms 


•Continuous surveillance for weather now - and 

FORECASTING/ DESCRIBING BACKGROUND KINEMATIC AND 
thermodynamic STRUCTURE OF THE ATMOSPHERE AND 
DESCRIBING AREAS IN WHICH DEEP/ MOIST CONVECTION 
IS INITIATED 


•Intense surveillance for inspection of atmospheric 

CHARACTERISTICS WITHIN AREAS OF OCCURRENCE OF DEEP/ 
MOIST CONVECTION 


•Measurements concerned mostly with broad scale 

PHENOMENA WITH EMBEDDED SMALLER SCALE PROCESSES 
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Early storm: 


•Concerned with small clouds and early phases of mature 

STORMS 

•Objectives 

SPECIFICATION OF ENTRAINMENT OR MIXING PROCESSES 

DETERMINATION OF PRIMARY ICE FORMATION REGIONS 

DEFINITION OF PRECIPITATION GROWTH TRAJECTORIES 

ELUCIDATION OF CLOUD EFFECTS ON AEROSOL AND CLOUD 
CONDENSATION NUCLEI POPULATIONS 

•Smaller spatial scales ('^100 's of meters) 


•Shorter time scales ('^1-2 minutes) 

•Multiple aircraft penetrations and glider ascent 
•Kadar detection of small particles and chaff for 

WIND DETERMINATION 
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VERTICAL SECTION 




^// Photographic 
Aircraft Track 
at Cloud Base 


^2* 10 km 
K ^ 



PLAN VIEW 



riATURE storm: 


•Characterized by tops in excess of '^9 km MSL or 

REFLECTIVITIES GREATER THAN '^^5 dBZ 
•Objectives 

. SPECIFICATION OF HYDROMETEOR EVOLUTION 

DESCRIPTION OF STORM STRUCTURE AND EVOLUTION 

DETERMINATION OF WATER BUDGET COMPONENTS 

AND FACTORS AFFECTING PRECIPITATION EFFICIENCY 

INVESTIGATION OF ENTRAINMENT OR MIXING PROCESSES 

•Spatial scales from cells to storm C'^2-50 km) 


•Temporal scales of '^2 min to storm lifetimes 
•Aircraft measurements in near environment and T-28 

PENETRATIONS 


•Radar mapping of precipitation and kinematic structure 
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Data management and analysis: 


•Field quality control 

•Level 1 data (from developed instrumentation) 

Reduced and archived to specified degree 
•Available to all participants as soon as possible 


•Level 2 data (from experimental instrumentation) 

Special reduction and limited archiving 
Availability negotiable with principal investigator 


•Case study approach to describe storm 


•Specific analyses for hypotheses testing 


E-24 


Operations and schedules: 


•Experimental period from early May through early 
AUGUST/ 1981 


•Overall coordination from operations center at 
. Miles City 

OPERATIONS director 
CP-2 SCAN controller 
Doppler radar coordinator 

AIRCRAFT COORDINATORS 

T-28 

ABOVE CLOUD BASE 
BELOW CLOUD BASE 


•Time table 

CCOPE OPERATIONS plan: Nov. 1980-Feb. 1981 
Participants meeting: Feb. 1981 
Distribution of operations plan: Mar. 1981 
Installation and calibration of observing 
systems: Mar. -Apr. 1981 
Rehearsal: MaY/ 1981 
Field program: May-Aug. 1981 
Removal of facilities and instruments: Aug. -Oct. 
. 1981 
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PART II OF PRESENTATION BY 
DR. GEORGE H, FICHTL, MSFC 
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A DATA PROCESSING 
V FLIGHT OPERATIONS 


APPENDIX G 


REPORT OF CONVECTIVE PHENOMENA TEAM 


The First Scientific Working Group Meeting 


MSFC Airborne Doppler Lidar V7ind Velocity Program 
NASA/MSFC, 25-26 August 1980 

Report of Convective Ph eno mena Team . 

H. Orville, R. Koenig, J. Miller, J. Telford. 

B. Jones, G. Alger, R. Lee, B. Bowdle. 

I. Introduction and Flight Plans 

A Group meeting was assembled to focus on the planning of specific 
experiments, to establish some priorities, identify interested 
scientists who would like to participate, establish any special 
requirements, make recommendations on data processing, and to prepare 
flight plan outlines. 

The types of experiment of interest to the group were discussed 
first, without regard to priorities or importance. After the ideas had 
been explored, it became clear that since the number of convective storms 
in the CCOPE (Cooperative Convective Precipitation Experiment) field 
experiment area would be limited to only a few days during the 
operational time period the flight plans had to be designed with a 
hierachy of abort experiments so that the easily identified and lowest 
probability events should take priority until their quota is filled. It 
was estimated from the basis of 60 flying hours total, from which 
testing and ferry time must be deducted, that optimistically 45 hours 
flying would be available in the CCOPE experimental area. With five 
hours per flight this would be 9 flights, frcxn which 3 would likely to 
be aerosol experiments which could not combine with the Convective 
requirements. This leaves 6 flights of 5 hours each for planning 
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purposes of the Convective team. 

The CV990 aircraft would probably be based at Ellsworth A,F,B. near 
Rapid City (210 miles to Miles City) or Minot A.F.B. (260 miles to Miles 
City) near Minot. And since it is desirable to cruise at 25000 ft. the 
pilots should have a good view of the cloud top situations from 10 
minutes or more before reaching the measurement area. Thus a flight plan 
to explore the outflow around Cumulo-nimbus anvils can be chosen at this 
time vhen such cloud conditions are present. If Cumulo-nimbus anvils are 
not visible at this stage, but we see clearly defined, growing, and 
decaying, cumulus turrets, below, say, 30,000 ft. then the entrainment 
study can be selected. 

Since it takes a 2 hour period to prepare for a flight, and 
probably an hour from the decision to take-off to reaching the vicinity 
of the experimental area, there may be occasions v*ien the cloud field we 
encounter on the site will be unsuitable for either of these measurement 
programs. In these cases a decision will have to be made as to whether 
to implement some study in the boundary layer, or to sacrifice the ferry 
time to site and abort the mission by returning to base. Another 
alternative is to seek targets of opportioiity beyond the CCOPE support 
area, and an\^ow it is possible that interesting measurements can be 
obtained during the ferry run to the target area. For example it is of 
interest to look at particle reflectivities in cirrus cloud. 

Flight patterns below cloud cannot be discussed in detail at this 
present time since tight cooperation with CCOPE operations will be 
essential for work in this area below cloud but we request such support. 
The study of storms gust front dynamics, another high priority mission 
we also wish to recommend, will be particularly sensitive to careful 
coordination. 
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The other experiment below cloud, to study the feeder flow for 
anall, non-precipitating cunulus, may be best carried out away from the 
storm of prime interest to CCOPE. In this case the aircraft will be 
below cloud so that the selection of suitable individual clouds will not 
be possible and long straight runs are probably the preferred 
observation mode. 

A very interesting data set could be gathered by flying encircling 
squares around a moderate sized growing cumulus at levels from above 
cloud tops to below cloud base. However since this lidar device requires 
clear air for its sixrcessful use there is some doubt as to whether there 
will be adequate useful opportunities, since the requirement is that a 
cloud is sufficiently isolated, from a base say at 10,000 ft., to, say, 
16,000 ft., so as to leave about 10 miles cloud- free air around it in 
all directions. If this occurs it should be taken as a target of 
opportunity. 

Other targets of opportunity should be accepted during instruments 
check out flights. If marine stratus is encountered near Moffat, or the 
clear air boundary layer, or cirrus, these would also be possible 
targets at this stage. 

2. Detailed Instrunentation Considerations and larget Area Flights 


Plans 

(a) Rectangles around clouds. 

The measurements will be taken in non-turning level flights except 
for the cloud- top entertainment studies. This requirement is set by the 
need to intersect every air parcel twice with a horizontal beam. If the 
two horizontal components of velocity are to be derived. It should also 
be noted that the air velocity is being measured relative to the 
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velocity of the instrument within the aircraft. Thus any relative 
movement between the optics of the lidar and the velocity reference 
given by the initial platform, will appear in the final result as an 
error, unless this displacement, and relative velocities and angular 
rates, are included in the processing operations to derive the true 
velocity. Rapid vibration is virtually impossible to account for in this 
way, and hence must be eliminated mechanically. Flexure of the airframe 
and mounts needs to be examined to be sure they are negligable at the 
hoped for accuracy, which we understand can be a fraction of a meter per 
second. The separation between the lidar and the platform will produce a 
velocity difference when the attitude of the aircraft is changing, 
heading is probably the largest component, and this must be corrected 
for. Thus inertial velocity (3 axis) and attitude (3 axis) must be 
recorded fast enough to allow smooth differencing to give rates covering 
the maximum response frequency of the aircraft, say, 2 cycles per 
second. Ihe vibration modes of the platform should be checked to ensure 
that no appreciable motion occurs beyond the recording frequency, since 
then aliasing would inject false motion into the data. 

A problem exists in obtaining precise verification of the 
instrunental performance since detailed matching of the air velocity of 
parcels in turbulent motion with, say, a radar chaff echo, is likely to 
be impossible, so that only smoothed averages can be checked. Thus if 
the atmospheric scientists using the data are going to have confidence 
in the measurements every precaution will need to be taken to 
demonstrate that we have accounted for all errors. No new discovery can 
be made with measurements which cannot be trusted. 

Each encircling flight to measure air motion around a cloud anvil 
will need to have flight legs parallel and perpendicular to the shear. 
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Since the lidar device points out from the left hand side of the 
aircraft only the velocities on the left hand side can be measured. 
Because turning-fl ight prevents the normal sampling density fore and aft 
it was the opinion of the convective team that all turns during sampling 
rectangles should be right hand turns of 270 degrees of arc, so that the 
straight and level flight would extend to the extremes of the area under 
study, which would often encompass a storm cloud and hence contain 
inaccessible areas because of cloud obscuration. With an estimated 
reliable range of 10 km (10 miles max.) it was felt that the opposite 
legs of the rectangle should be about 20 km apart, provided the cloud 
fitted comfortably within this dimension. In enclosing a cirrus anvil 
the legs heading down shear would usually need to be extended, perhaps 
to 100 km on occasions, to examine a large part of the ice crystal 
region. 

There is a real question as to whether the scattering aerosol will 
be dense enough at the higher levels to give lidar returns. If the 
returns are too uncertain these data gathering ideas will have to be 
revised. 

For the rectangles below the level of the anvil it is important to 
examine the air in the direction down shear, but below the anvil, where 
the air is free of ice crystals. This air will probably show signs in 
the moisture and velocity fields of having been part of the cloud a 
little earlier. Thus these rectangles should cross beneath the anvil as 
close to the cloud as the pilots feel can be accomplished without risks 
of encountering hail. Modifications of these plans may have to be made 
on site to avoid precipitation, or flying too much in surrounding cloud, 
vhich will obscure the lidar. 

It is a requirement for the software to clearly distinguish the 
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cloud edge in a lidar return and ensure that the cloud edge is clearly 
delineated on all plan plots, such as that of velocity. It was felt that 
the proposed velocity vector field was too complex to be assimilated 
\4iere the tag on the wind vector needed to be interpreted for both 
length and angle. It was suggested that a better display would be to 
have a simple arrow for wind magnitude and direction with a blob as an 
arrovJiead, where the area of these elements represented range corrected 
echo intensity (and to display velocity variance separately) . The aim 
would be to have the density of an area proportional to its 
reflectivity. 

It was generally felt that an experimenter needed to receive a raw 
data tape with all the data merged to a single tape which he could 
readily read into a Fortran (or other) language program. 

(b) Cloud top Study 

The second level experiment to be chosen in flight is based on the 
need to study motion of the cloud surface at the top and sides of cloud 
turrets where it is thought the entrainment of dry air into the cloud 
begins. Since, ideally, we need the vertical component of the motion of 
the cloud edge, because the buoyant energy release transfers initially 
to vertical motion, a single doppler lidar beam looking downwards would 
be the first choice. However since there are theoretical reasons for 
expecting that the component of velocity responsible for the entrainment 
is the component perpendicular to the cloud surface, and since the 
normal anall wandering of the aircraft heading will prevent locating the 
beam on any particular spot of cloud, there seems no reason why a beam 
pointed down at 45 degrees of arc below the horizon while the aircraft 
circles above cloud top should not give extremely valuable information. 
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To begin with, until the operation and cloud behaviour are better 
understood, this is probably the preferred choice, anyhow. 

Thus the flight plan calls for the aircraft to circle above the 
cloud turret. Since the lidar beam can be deflected downwards by 20 
degrees by the optics, then 30 degree banks on the aircraft would be a 
suitable flight parameter for determining the height and turning radius. 
Some sort of crude optical siting device in the cockpit could be used to 
assure the pilot that he was maintaining radar contact with the chosen 
cloud turret top. 

In this operation we are seeking to obtain a detailed knowledge of 
how the surface of the cloud responds to the in-cloud turbulence, and 
how this results in exterior dry air entering the cloud through this 
surface boundary. Thus we would request that the capacity of the data 
system be devoted to recording as much information as possible about the 
regions of penetration just into the cloud surface, and from the air 
just outside. The fore and aft scanning procedure would need to be 
stopped, and by using the computer to range gate over say 10 to 20 bins 
at the cloud edge, atxl using individual 1 micro-second pulse gates 
without any averaging, the full data capacity could be devoted to this 
narrow range of interest. This range selection would need to be a 
continuous software monitored function. It is also desirable to use 
software amplitude control to avoid saturation of the cloud edge, so 
that diluted cloud regions, where evaporation has reduced the liquid 
water content, can be faithfully recorded for intensity of echo as well 
as for their velocity profile. 

It would be nice to scan across the cloud surface so that 
successive radar hits were contiguous, but this is too much to expect, 
except by relying on chance scans like this every now and again in the 
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records. It is essential to have a video camera coaxial with the lidar 
beam to record the cloud surface visually in real time as these flights 
proceed. Visable time on each frame is needed. In this way the wander 
velocity of the cloud image should be visibly evident. With care it may 
be possible for the pilot to maintain the cloud image rotating smoothly 
in the beam, so we measure an advancing track around the cloud turret. 

It is not clear what duration of this circling flight can be 
maintained with accuracy and without incapacitating passengers, we would 
suggest one hour for planning purposes. A total period of four or five 
hours in the season is a reasonable target. 

The subcloud, non-precipitating cumulus, feeder flow study 
similarly needs special operational considerations. The problem here is 
that when the aircraft is a few hundred feet below cloud base in an 
extended area of cloud it is not possible to see the cloud tops or 
sides. Thus the desired flight profile of a rectangle located so as to 
scan the area below a small vigorously growing cumulus is probably 
impractical since we cannot select the position of such clouds when we 
are beneath them. Thus the practical approach is to fly long legs 
beneath the cloud region, up and down the shear direction as based on 
the subcloud and mid cloud wind difference, and record on the video tape 
all the information that can be seen of the cloud bases. In this way 
cases can be selected where a good correlation between the wind field 
and its location below cloud base will let us establish the flow 
structure relative to the cloud. Ihis needs to be done at several levels 
below cloud base on each occasion, say 500 ft., 1,000 ft. and 1,500 ft. 
below bases. The runs would then be perhaps 30 to 60 km in length with 
turns at the end to bring the return back 10 km. to the left so that the 
same area is scanned again on return. The altitude for these runs will 
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probably be 10,000 to 16,000 ft. and will require careful coordination 
with aircraft control. 

The gust front study may be possible as an airborne alternative 
choice, but probably will need separate selection as a mission because 
of the tight coordination needed for the aircraft. The plan would be to 
descend beneath the storm or squall line and locate the front with 
guidance from the radars or other aircraft. It is then desired to fly 
parallel to the front a few kilometers ahead of it, with a return track 
on roughly the same line as the previous track so as to complete the 
circuit and begin the next leg. The legs would typically be about 30 km 
for a single storm, and 100 - 500 km for a squall line, in length, with 
the procedure turns beyond this. Depending on terrain it would be 
desirable to repeat the pattern at altitudes of 100 m, 300 m, 1,000 m, 
and 2,000 m, or as time permits. In rough terrain these heights would 
need to be modified. 

If all weather failed, then the ultimate backup experiment would be 
to study the aerosol contents and velocity distributions in the 
convective planetary boundary layer. This would simply involve a race 
track box pattern at various levels from just above the inversion so 
near the surface. A similar procedure at higher altitude would enable 
cirrus particle densities and motion to be studied with altitudes from 
just above cirrus tops to just below the cirrus base. These studies 
could also be condixrted on long iiorizontal stretches during ferry runs. 

The marine or continental stratus observations would be similar to 
the bouiidary layer aiod cirrus layers mentioned in the last paragraph 
except that the cloud top banked circle (anticlockwise) entrainment 
observations should be added, and occupy a substantial time period, 
since these are the most important observations in this case, and 
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flights at levels within cloud are not likely to be useful. Subcloud 
turning observations (clockwise) below cloud base to measure vertical 
motion of the base of marine stratus would be extremely valuable when 
adequate airspace is available. 
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THE FIRST SCIENTIFIC WORKING GROUP MEETING 
MSFC AIRBORNE DOPPLER LIDAR WIND VELOCITY PROGRAM 
NASA/MSFC, August 25-26, 1980 


Report of Secondary Flows, Boundary Layers, Turbulence , 

and Wave Team 

Team Members: J. Scoggins, D. Fitzjarrald, R, Doviak, 

W. Cliff, (W. Frost absent) 

1 . Introduction 

The following is a brief summary of the report to the entire group 
by the boundary layer team. The primary objects of discussion were pos- 
sible test plans and their relative merits and operational difficulties. 

Four likely candidates for the flight tests were discussed in some detail 
and are given later. General criteria for a flight test option are that:. 

(1) there be a good opportunity for comparison with other measurement 
techniques, (2) the flow to be measured is of considerable scientific or 
practical interest, and (3) the airborne laser Doppler system is well 
suited to measure the required quantities. The requirement for comparison, 
i.e., ground "truth", is particularly important because this will be the 
first year of operation for the system. It will be necessary to demonstrate 
that the system does actually measure the winds and compare the results with 
other methods to provide a check on the system error analysis. The unique- 
ness of the laser-Doppl er system precludes any direct comparison, but 
point measurements from tower-mounted wind sensors and two-dimensional fields 
obtained from radars with substantially different sampling volumes will 
be quite useful . 

The flight test options presented below can be economically grouped 
to minimize ferry time. The first two (California Central Valley and 
San Gorgonio Pass) are in California, so that the aircraft can be based 
at NASA Ames. They should be done early in the program, as soon as the engi- 
neering tests are completed. The next test would logically be at NSSL in 
Oklahoma, where the schedule is flexible and opportunities for tower and 
radar comparison are good. By this time considerable operation experience 
should have been gained, so that moving up to Montana to the CCOPE experiment 
would be the next step. The CCOPE is by far the most complex operationally, 
requiring integration into a very big and busy experiment. Success would 
be maximized by having as much experience as possible with the operation of 
the system and some confidence that the system does measure what it is 
supposed to measure. On the way east to NSSL and back west from CCOPE, the 
Boulder NOAA/NCAR tower instruments can be compared in a fly-by. 

In addition to the flight test options, the real-time data display 
was discussed in some detail. The recommendations are presented later. 

A summary of the flight test options is given first, followed by the dis- 
cussion of the data display and, finally, by details of the flight test 
options. 
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Summary of Flight Test Options 
A. California Central Valley 


Purpose : Investigate the spatial and temporal variation of the 

boundary layer wind flow patterns in the Central Valley, examine the detailed 
flow patterns in the vicinity of ag-burning smoke plumes, and the local flow 
near the Geysers geothermal field. 

Data Comparison : 1500 ft tower at Walnut Grove, possibility of tracer 

release by DOE Livermore, surface wind field at the Geysers. 

Flight Requirements : Racetrack pattern around Central Valley near top 

of mixed layer (approximately 3000 ft), continue from noon through evening. 
Smoke plumes as targets of opportunity, box patterns at various levels at 
the Geysers. 

Interest : Important problem in applied meteorology, of great interest 

to air pollution researchers. Interesting and complicated local flow, 
very little data available. Collaboration from California Air Resources 
Board, EPA, DOE Livermore, and DC Davis. 

Suitability : Well suited for measurement by Doppler system, with 

space scales and resolution within system possibilities. Good place to try 
out the different resolutions, etc., of the system. 

B. Wind Resource Assessment of San Gorgonio Pass, California 

Purpose : Investigate spatial and temporal variations of flow through 

mountain pass in Southern California. 

Data Comparison : 100 m tower, 50 m tower, twelve 10 m towers, and 

several years of averaged data. 

Flight Requirements : Back and forth pattern at several levels at 

the end of the pass. Continue for sufficient time (approximately 3 or 4 hours) 
to see the evolution of the flow. 

Interest : Prime candidate for wind energy farm. Spatial details of 

the flow needed to compare with surface and small tower data. Climate type 
wind data available for 3 years, but nothing to give the extent in space of 
the wind resource. Collaboration with Southern California Edison, DOE 
Wind Energy Program, California State Wind Energy Commission, Battelle Pacific 
Northwest Laboratories. 

Suitability : Space and time scales well within system capabilities. 

Interest is in the fine-scale details of the motion (approximately 0.5 km). 


H-2 


C. Radar/Tower Comparison at NSSL 


Purpose : Fly-by 1500 ft tower in Oklahoma City and cover the same 

area as NSSL dual-Doppler radars. Compare the measurements by two dif- 
ferent instruments and determine the two-dimensional spectra of the wind 
at different levels in the boundary layer. 

Flight Requirements : Fly-by tall tower. Racetrack pattern at different 

heights within the boundary layer. Pattern coinciding with area covered 
by NSSL radars. 

Interest : Excellent opportunity for comparison with another instru- 

ment that measures a two-dimensional wind field in clear air. Wind 
patterns and spectra can be compared. Scientific interest in horizontal 
spectra. Collaboration with NSSL. 

Suitabitity : Space and time scales within system possibilities. 

Interested in all the scales within the flight box that can be resolved 
by the system; therefore, the finest possible resolution would be desired. 

D. CCOPE, Large Field Experiment in Eastern Montana 

Purpose : Compare boundary layer wind measurements with those meas- 

ured with radars (clear air and chaff release), aircraft, and surface 
measurements at the CCOPE test site. Measure the before-cloud and under- 
cloud convergence in the boundary layer. 

Flight Requirements : To be coordinated with CCOPE to avoid conflicts 

with multitude of other aircraft operating in area. Will schedule some 
flights on "off days" or early in the day to be able to compare with radars 
without getting in the way of other aircraft. 

Interest : Opportunity for contribution to large experiment. Many 

collaborators. Should be confident of the system and experienced in opera- 
tion because this is not the best place to experiment with it due to large 
number of other aircraft. 


Suitabil ity : Same as for boundary-layer measurements at NSSL. 

3. Real-Time Data Display 


As a result of discussion within the group and with the M&S personnel 
who designed the display software, it was concluded that the intensity (S/N) 
information was not adequate. It was not possible to look at the vector 
display and see the intensity information easily. It was believed that this 
would be a most important quantity to look at during the flight. There was 
considerable discussion regarding the form of the display, but no clear 
consensus was reached. Any display oT Uie scalar field of intensity would 
be easier to read than the arrowhead lengths on the wind vector display. 


4. Details of Flight Test Options 

Some details of the proposed flight tests are given below. In two 
cases they represent a summary of the group discussions. For the San Gorgonio 
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Pass test the description is supplied by W. C. Cliff. Detailed flight plans 
were discussed by the group, but in view of the many uncertainties it appears 
better to proceed in a more general manner at this point. Certainly, it will 
be desired to try all the different resolutions and operation modes in the 
first tests. Operation at different turbidity levels will also be necessary 
to see what range can actually be attained. The first tests should provide 
all these opportunities, together with ground truth for comparison. Assum- 
ing that all is in order, the last test should be to contribute to the CCOPE 
experiment. 

California Central Valley Flight Test Option 

Justification ; The air flow and circulation patterns within the 
Central Valley are of great interest because of their importance in air 
pollution assessment and control. Pollution within the valley at present 
is primarily due to agricultural burning, but the increasing pressures of 
urbanization and' the proposed siting of a large fossil-fuel power plant 
indicate that other sources will dominate in the near future. The wind 
data available at present are not adequate to assist the regulators who are 
controlling agricultural burning (California Air Resources Board, various 
county air quality groups) or the planners interested in long-term prob- 
lems (EPA, ARB, electric utility, DOE). A high-qual ity measurement of the 
temporal and spatial variations of the flow would be a significant contri- 
bution. 

In addition to the general flow patterns, there are two additional 
flow measurements that would be of considerable use. By studying the details 
of flow near the large ag-burning smoke plumes, it may be possible to mini- 
mize the pollution. Some work on this has been done by researchers at 
UC Davis, but additional data would be helpful. A large interagency effort 
has been conducted to determine the details and climatology of the local 
air flow in the vicinity of the Geysers geothermal field in northern 
California. A short-term measurement of the wind flow patterns would be 
useful to the mesoscale modellers (UC Davis, DOE Livermore). 

This area provides a number of advantages for the testing of the 
airborne laser-Doppler system. It is close to NASA Ames, has some possi- 
bilities for comparison of data (1500 ft tower at Walnut Grove, instrumen- 
tation at the Geysers, tracer release by DOE Livermore), has at least one 
place (Sacramento Delta) where the flow is quite regular and known in direc- 
tion, and is normally quite high in optical scatterers. 

San Gorgonio Pass Test Option 

Principal Investigator : Dr. William C. Cliff, Battel le PNL 

Objective ; To characterize the vertical and horizontal extent of the 
accelerated flow region through and on the east side of San Gorgonio Pass. 
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Location ; San Gorgonio Pass, California. San Gorgonio Pass is located 
approximately 120 miles due east of Los Angeles, California. 

Description of Area ; San Gorgonio Pass is approximately 40 km (25 miles) 
long with a westerly elevation of approximately 760 m (2500 ft) gradually 
dropping to an elevation of approximately 200 m (700 ft). The width of the 
pass is generally approximately 8 km (5 miles) [at a contour approximately 
300 m (1000 ft) above the pass floor]. On the east end the pass quickly 
broadens and becomes an open desert floor. The mountain rises to approxi- 
mately 3350 km (11,000 ft) within 16 to 24 km (10 - 15 miles) on each side 
of the valley floor. 

Regions for Wind Assessment : The easterly 8 km (5 miles) of the pass 
and the adjoining 16 km (10 miles) of desert valley floor. The valley floor 
area to be assessed is roughly 16 km (10 miles) in the east-west direction 
by 24 km (15 miles) in the north-south direction. The pass area to be 
assessed is roughly 8 km in the east-west direction and 8 km in the north- 
south direction. 

Elevations of Horizontal Wind Field Mapping ; Anticipated flight paths 
should be at approximately 100 m (0.330 ft), 300 m (MOOO ft), 500 m (M600 ft), 
and 1500 m (<^^5000 ft) above grade level over the valley and with the potential 
addition of a 1000 m ('^^3000 ft) elevation taken in the pass (refer to Fig. 1). 

Desired Horizontal Spatial Resolution of Velocity Field : For ele- 

vations below 30o m (approximately 1000 ft) a resolution of < 300 mis 
desired. Above 300 m a resolution of • 500 m is desired. 


Ground Truth Data Availabil ity: At the expected test period the fol- 
low i ng~meteoroTogTcaT~data^ will be operating and the data will be 

available for the test program; 

a. One 100 m (330 ft) tower, 4 levels of instrumentation 

b. One 50 m (160 ft) tower, 3 levels of instrumentation 

c. Approximately twelve 10 m (33 ft) towers with instruments at 
10 m level only 

Previous Data Available : 

a. One 100 m (330 ft) tower; approximately 1 year of data 

b. One 50 m (150 ft) tower; approximately 3 years of data 

c. 20 stations at 10 m (33 ft); 1 year (August '78 - August '79) 

d. 12 stations at 10 m (33 ft); approximately 1.5 years 
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Agencies Expressing Interest in San Gorgonio Test Option ; 

a. DOE Federal Wind Energy Program/Wind Characteristics Program 
Element 

b. Southern California Edison Co. 

c. State of California/California Wind Energy Commission 

d. Battelle, Pacific Northwest Laboratories 

Justification ; The San Gorgonio Pass region is a known high wind 
resource region and is currently identified as a DOE candidate wind turbine 
site for the potential placement of one or multiple large wind turbine(s). 

The area is of extreme interest to the State of California Energy Commission 
and the Southern California Edison utility company as a high wind resource 
region for the potential placement of a large wind turbine farm for electri- 
cal power production. Southern California Edison has already purchased 
two large wind turbines from private companies, the first of which has been 
field erected and is expected to become operational later this year (1980). 

To establish the number of large wind turbines that this region's 
wind resource could support, the vertical and horizontal extent of the 
accelerated flow resource needs to be determined. The San Gorgonio Pass 
accelerates the flow from the west side of the mountains through the pass 
and spreads the accelerated flow onto the desert floor on the easterly side 
of the pass. The vertical extent of the accelerated flow is required to assess 
the potential mountain flow which may feed to the lower levels if wind 
turbines extract energy from the lower level winds. The vertical extent is 
needed to determine how much land comprises the rich wind resource area for 
wind turbine placement. There is extensive ground truth data assured during 
the testing, as described under "Ground Truth Data Availability," for system 
verification of the airborne NASA Doppler lidar system. This program pro- 
vides a direct spin-off of advanced NASA technology to other government 
agencies, DOE, and the State Energy Office of California, a private utility. 
Southern California Edison, and to the general public who will benefit from 
the power developed by optimally placed wind turbines in the test option area. 

NSSL, Oklahoma Flight Test Optio n 

Justification ; This flight test provides an excellent opportunity 
to test the airborne laser Doppler system against another type of system that 
can measure a two-dimensional field of wind in clear air. NSSL has already 
used their dual-Doppler radars to look at the convection in the afternoon 
boundary layer. The laser system will provide a finer scale result to be 
compared with the radar-measured winds. The laser resolution is approxi- 
mately a factor of ten smaller. The horizontal spectra can be compared in 
the larger scales, and the laser data will provide spectra of smaller scale 
turbulence. A well-instrumented 1500 ft tower is also available to check 
one point on the laser -measured wind field. Data reduction from the radars 
can be done quickly, so that a nearly day-by-day check can be made that the 
two systems are in agreement. 
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A significant advantage of this test option is that there is no schedul- 
ing problem. The radars are installed at NSSL and can be operated whenever 
the 990 is ready to fly. The tower instruments are also in regular operation. 
NSSL has conducted aircraft operations in both areas before, so there should 
be no difficulty fron an operations standpoint. The absence of a large field 
program is an asset in this because it will allow the 990 laser system to 
set the schedule. This should be ideal for the first-year test, in view of 
anticipated difficulties in operating a new system and gaining experience in 
coordinating the aircraft movements to obtain the best possible wind data. 

Participation in CCOPE 

The CCOPE project is a very large field experiment aimed at cumulus 
dynamics on the high plains (eastern Montana). A large number of aircraft, 
Doppler radars, ground-level winds, and radiosondes are involved. Planning 
is elaborate, and coordination will be difficult. The project offers a good 
opportunity to make a contribution with the airborne laser-Doppler system. 

In the boundary layer, convergence and wind patterns in the clear air before 
storms develop and underneath the developing storms will be of great interest. 
At upper levels, convergence of aerosols and cloud-edge motions will be of 
interest. Such a complicated experiment is not the place to try out a new 
system, however., The airborne laser Doppler should be used enough before 
starting this project to ensure that the system is actually measuring what 
we believe it to be measuring. The operating procedures (flight tracks, 
altitudes, spatial resolution, range of operation, etc.) that yield good 
wind fields should be well proved by the time of the CCOPE. This will ensure 
that the participation will have a high yield and that the cooperation from 
the CCOPE participants will be worth their effort. 

There will exist opportunities for testing the laser system using the 
radars that will be in position in the same manner as at NSSL. Either in 
clear air (a few of the radar sets) or with chaff drops (most of the radars), 
two-dimensional wind fields can be measured by the radars and compared with 
the laser system. Such comparisons will have to be on "off days" of the 
big experiments or early on flight days. This would entail some problems 
for the CV-990 because of the time necessary to plan flights and the expense 
of waiting in readiness for an appropriate time to fly. 

Because of the scheduling difficulties inherent in such a big experi- 
ment, it appears a better idea to accomplish most of the operational and 
comparison testing in prior tests (the ones in California and at NSSL offer 
good opportunities) and to go to the large experiment with a system that 
is proven. 
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Report by Team Member 

Dr. Richard Doviak, NSSL, Norman, Oklahoma 


U.S. DEPARTMENT OF COMMERCE 

National Oceanic and Atmospheric Administration 

ENVIRONMENTAL RESEARCH LABORATORIES 

National Severe Storms Laboratory 
1313 Halley Circle 
Norman, OK 73069 


September 9, 1980 


Mr. George Fichtl, ES82 

George C. Marshall Space Flight Center 

Huntsville, AL 35812 

Dear George: 



Thank you for inviting me to the first scientific working group meeting of 
an airborne lidar wind measurement program. We at NSSL look forward to a Joint 
NASA-NSSL experiment utilizing NASA's airborne Doppler lidar and NSSL's Doppler 
radars. I take this opportunity to confirm that our Doppler radar should be 
available for the joint experiment late June or July 1981 and we are attracted 
by the prospect of having the opportunity to compare the kinematic structure of 
the atmosphere measured independently by these two systems. 


Upon my return to Oklahoma, we thought of some other experiments which you 
might want to consider conducting during the time the CV 990 is in Central 
Oklahoma. Hopefully, the aircraft can be stationed at Tinker Air Force Base for 
a week to 10 days during which time we could complete the listed experiments. 
They are: 


I. Comparison of Horizontal Wind Fields 

II. Comparing Two Dimensional Spectra of the Kinematic Structure of the 
Convectively Driven Planetary Boundary Layer 

III. Comparison of Doppler Derived Heat Flux Profiles with Those In Situ 
Aircraft Measurements 

IV. Comparison of Mesoscale Divergence and Deformation Measurements in the 
Pre-Storm Environment 

V. Comparison of Thunderstorm Gust Front Winds 


I . Comparison of Hori z ontal Wi nd I- i el ds 


This is one of the experiments we discussed at the meeting and is one that 
should be relatively easy to execute. The aircraft could fly at a fixed altitude 
(e.g., 1 km) some 40 km SW of Norman, and the radar and lidar both should have 
sufficient tracers on any sunny afternoon to map the wind. In this experiment 
we will compare (1) the mean horizontal wind averaged along the strip (e.g., 

10x30 km*" of atmosphere probed by Lhe lidar, and (2) the secondary wind field 
patterns. I'm hopeful that we will get at least good agreement for (1) and 
comparisons could be made shortly (1-2 days) after the experiment. If we get 
good agreement with (2), I'll be delighted. In this experiment the aircraft 
should fly crosswind and along wind and comparisons of mean and secondary fields 
made. 
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II. Comparing Two Dimensional Spectra 

, Thi S' experiment would provide data for an analysis of the type described in 
the enclosed paper, "Turbulence and Waves in the Optically Clear Planetary 
Boundary Layer Resolved by Dual Doppler-Radars" . Of particular interest here 
wQuldpbe to extend these measurements to obtain vertical profiles of u, v, w, 

^u’ '^v’ >"adar and lidar. This experiment could be conducted as part of 

Experiment I, but it would be necessary for the aircraft to make passes at 5 
different altitudes (e.g., 300, 600, 900, 1200, and 1500 m). Can the aircraft 
make direct measurement of w and This would be important because computation 
of w from lidar and radar data is much more difficult, and it would be imperative 
to have a direct in situ measurement of turbulence even if it is along a single 
line whereas the lidar and radar provide an area of data. 

III. Comparison of Heat Flux 

I've enclosed a copy of a paper, "Measuring Heat Flux and Structure Functions 
of Temperature Fluctuations with an Acoustic Doppler Sodar"★ ** There it describes 
a technique that can be used to relate profiles of vertical wind variance to 
profiles of heat flux. Although we have not performed the radar measurement 
yet. I'm confident that with minor modifications to our radar, we should be able 
to obtain the requisite data. We have a shrouded vertically pointed antenna 
that makes sidelobe levels small and I believe we would be able to collect data 
as low as 300 m above the surface. In this experiment we hope that the aircraft 
would be able to make spiraling ascents around the vertically pointed Doppler- 
radar if measurements of w are not 2 Compromised significantly. What we need to 
compare first is the profiles of and if the aircraft is suitably instrumented, 
the heat flux profiles. Can the aircraft measure directly heat flux (by correlation 
of fluctuations in w and 6)? As an adjunct to this experiment, we may want to 
perform fly-by's over our meteorologically instrumented 444 m tower where we can 
make heat flux comparisons with tower measurements. 

One of our researchers (Bob Rabin) has recently suggested (see 2 attached 
memo for the record) that the refractive index structure constant C^ can be 
obtained from knowledge of sensible heat H and latent heat E fluxes. Can the 
aircraft provide data so we can check this hypothesis? Can one make water vapor 
variance and temperature covariance measurements with the aircraft's radiometer? 
What are the characteristics of the radiometer? 

IV. Comparison of Mesoscale Divergence and Deformation in the Pre-Storm 

Environment 


This may be the most trying experiment because we would have to wait until a 
storm situation develops in Central Oklahoma. Nevertheless, if such conditions 
do occur while the aircraft is here, I suggest we conduct it. I enclose a 
summary of a paper entitled, "Statistical Considerations in the Estimation of 
Wind Fields from Single Doppler Radar and an Application to Pre-Storm Boundary 


★ 

This paper is not included in this report because of copyright restrictions; 
please see Radio Science, 15, No. 2, March-Apri 1 1980, pp. 297-317. 

icic 

This paper is not included in this report; please see Journal of Applied 
Meteorology, 19, February 1980, pp. 199-205. 
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Layer Observations", that suggests that a single Doppler radar can map the 
mesoscale (L>10 km) patterns of divergence and deformation which we believe are 
important in the triggering of thunderstorms. This work is new at NSSL and 
untested. We suggest that the aircraft 2 fly triangular patterns at constant 
altitude over large areas (say 30x30 km^) and divergence and deformation computed 
using the mean winds measured along the path. 

V. Comparison of Thunderstorm Gust Front Winds 

This experiment will be the most difficult to perform because it requires 
the presence of a not too frequent phenomena within our radar range (=^60 km) and 
good coordination between aircraft and radar. The forecast of the gust front in 
the time (2 hrs.) it takes the aircraft to be airborne from a standby position 
makes difficult the success of this experiment. Nevertheless, wind shear measurements 
in the clear air outflow of thunderstorms is important in assessing the capabilities 
of Doppler radars to determine gust front hazards to landing aircraft. In this 
experiment we would direct the aircraft, using radar observations, to fly along 
the gust front so that wind in the gust front measured by radar and lidar can be 
compared. 

George, I hope I haven't overwhelmed you with experiments, but if we can do 
any one of these, I would consider the mission a success. Experiments I, II, and 
III could be accomplished in one day, IV and V in another. I would appreciate 
any information you would send me that describes the aircraft instrumentation 
and the products available. Give my regards to Ed Weaver and Jim Bilbro and 
please thank John Kaufman for kindly transporting me to the meeting. I didn't 
realize my motel was out of the way of his normal commuter route. 

Sincerely, 

Dick Doviak 

Chief, Advanced Techniques 


Ends. 

cc NSSL Managers 
Lee 
Zrnic 
Rabin 
Raj an, OU 
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U.S. DEPARTMEIUr OF COMMERCE 

National Oceanic and Atmospheric Administration 

ENVIRONMENTAL RESEARCH LABORATORIES 

National Severe Storms Laboratory 
1313 Halley Circle 
Norman, Oklahoma 73069 

August 27, 1980 


Memo for the Record 
FROM: Robert Rabin 

SUBJECT: Using radar to determine heating rate and evaporation near the 

earth's surface for an unstable boundary layer 



Purpose: We would like to ascertain as much information as we can about the 

prestorm environment using Doppler radar. Besides the obvious use of Doppler 
derived wind fields, radar reflectivity from refractive index fluctuations 
may provide information on the temperature and moisture structure of the 
boundary layer. Stephen Burk and Dusan Zrnic have already suggested how 
reflectivity could reveal a weakening of the moisture capping inversion 
(see Memo for the Record - Dusan Zrnic, 12 May 1980). The purpose of this 
note is to review what is known about reflectivity in clear air as related 
to the turbulent fluxes of momentum, heat and moisture in the planetary 
boundary layer. 

In clear air, the radar reflectivity, n (length”^), can be derived 
from the refractive index structure constant, Cjf (length "2/3) and radar 
wavelength, X (equation #1); see Hardy et al . , 1966, J. Geophys. Res., 71, 
pp 1537-1552. 

n - 0.38 (1) 

(Refractive index, n, is nondimensional ) 

Recall also from Gossard 1977, Radio Sci., Vol . 12, No. 1, pp 89-105, 
equation #2. 


.2 2 2 . .2 

Cp - a Cq + b 


s' - ^'eq 


where Cg^ (Deg 'length"^^^) , Cq^ (Water vapor density -length” and 
f-eq (Deg -water vapor density- length“2/3) are the structure constants for 
potential temperature variance, water vapor variance, and temperature - 
water vapor covariance, respectively, a and b depend on the mean air 
temperature, water vapor density, air pressure, and radar wavelength: 


- 30 

are typical : 


-1\ ._3n 

>• 


(m^gr" ). For 10 cm radars the following constants 


a^ = 2.24 X 10‘ 
b^ = 17.8 X 10" 
2ab = 12.6 X 10" 


Tropical Maritime Air 
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a^ = 1 .25 X 10"^^ 
= 16.2 X 10"^^ 
2ab = 9.01 X 10’^^ 


Warm Continental Air 


The radar detects backscattering from refractive index fluctuations of 
5 cm scale size. The associated temperature and moisture variations are in 
the inertial subrange. There, the structure constants are functions of 
molecular dissipation rates of turbulent kinetic energy e, temperature 
variance co, water vapor variance £ q , and temperature-water vapor covariance 
Egq (Equation 3); see Wyngarrd, et al., 1978, JAS, pp 47-58. 

C^e= 4 A (3a) 

cV 4 A e E (3b) 

= 4 

where these are nondimensional constants Ae = 0.4 

Aq = 0.4 
A0q = 0.5 - 0.6 


Combining equations 1, 2, 3, the radar reflectivity can be written in 
terms of the dissipation rates (equation 4),. 


n = 1.52X "''/V‘'/^(a^AQeg + b^A^e^ - 2 ab e^q) (4) 

The dissipation rates appear in the budgets (balance equations) of 
turbulent kinetic energy, temperature variance, water vapor variance, and 
temperature-specific humidity covariance. Fortunately, these budgets are 
best understood for unstable conditions (which concern the prestorm environ- 
ment). Observations of the individual terms in the budgets suggest that the 
dissipation rates are functions of the following surface layer* parameters: 

1) Momentum Flux: « « 

u^ = -u’ w^ = t/po (length /time ) 

where t ® shearing stress 

u*= friction velocity 

u'= fluctuation in wind component parallel to mean wind 

w'= fluctuation in vertical wind component 
overbar represents ensemble average 

Po= surface air density (mass/volume) 

2) Sensib le hea t flux: 

H=fbCpo'W (Energy/unit area/time) 


*The surface layer refers to the lower 10' s of meters of the planetary 
boundary layer where the fluxes are nearly constant with height. 
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where Cp = specific heat of air at constant pressure (Energy /mass/deg) 

e'= fluctuation in potential temperature (Deg) 

3) Latent h eat flux: 

E = La q'w' (Energy/area/time) 

where La = latent heat of vaporization (Energy/mass) 

q' = fluctuation in water vapor/density (mass/volume) 


The equations 5 a-d are based on surface layer measurements with tf 
limitations of horizontal homogeneity and stationarity. (See Champagne 
et al., JAS, 1977, pp 515-530.) 


e = 


KZ 


(1 


+ 0.5|^ 


2/3 3/2 

I ) 


where K = 0.35 (universal constant, nondimensional ) 


(5a) 


Z = height above ground (m); Z>Z 

2 ^ 

Zq = roughness length (0) - 10" m 


L = [-u^ Tv/Kg(H/pcp + 1.72 X 10"®(T)^i|)] 
Tv = Mean Virtual Temperature (“K) 

T = Mean Temperature (°K) 

O 

g = Gravitational Constant (m s" ) 

2 . 

KZ( pCp) ir; 


(5b) 


v/here 

u 


e 


q 


= nondimensional vertical gradient of potential temperature. 
= (1/1.35)* (l-of)"''/^ 

KZ(La)^u^ 


where (|) 


q 


e 


Oq 


nondimensional vertical gradient of water vapor density 

(l/l . 35 )' (1 - 9 ^ )V2 


2*H*E-»q 

KZ(pCpLa)u* 

(See Wyngaard, 1978, JAS pp 47-58) 


(5d) 
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Owing to experimental difficulties in measuring humidity fluctuations 
accurately, equations (5c-d) are perhaps less reliable than desired. However, 
all the equations are semi -empirical and should be considered approximations 
which best fit experimental data. 

The Doppler radars cannot usually measure reflectivity accurately so 
close to the ground to concern the surface layer. Hence, we must consider 
parameterization of the dissipation rates above the surface layer. Such 
parameterization is based primarily on AMTEX observations. Equations 6(a-d) 
follow from Lenschow and Uyngaard, JAS 1980, pp 1313-1326; Wyngaard et al . , 
1978, JAS, pp 47-58; and Lenschow, JAS, 1974, pp 465-474. 




0.57 

(SP+3.75) 


(SP-SP)+0.7(B-B)] 


(6a) 


where B = 1 - 1.5 1/1\ is the normalized buoyant production 

of turbulent kinetic energy (nondimensional ) 


- 

B = / BdZ = 0.25 
Z=0 

Zi = inversion height (m) 


SP = - [1 - 15^-]’^/'^ 


is the normalized shear production of 
turbulent kinetic energy 


SP 


= 21n + In i-- ^ - 2tan‘^X+ I dZ 


Z=0 


X = (1 - 15Zi/L) 


I/4 


S"W W h -2)+9.3(l-H ] 


^i 


(6b) 


where D= 3.059 Following Wyngaard et al . 1978 


when §= 0.61 ^ 


or D = 1 .4 


following Lenschow et al . 1980 
u 96 1 /I ^ H r 


3Z 
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'La' 'PCp 






.353 




] 


(6c) 




(6d) 


Equations 6 do not apply at the upper limit of the PBL (Z> .8 Zi ) 
because of entrainment through the inversion layer which is not only related 
to surface layer parameters. Other factors such as localized wind shear and 
wave phenomena also strongly affect the budgets of temperature and humidity 
variance and turbulent kinetic energy in that region. 

Results of higher order turbulence models give u* as an implicit 
function of surface layer geostrophic wind (pressure gradient) and sensible 
heat flux; equation 7 from Arya-Monthly Weather Review - Feb. 1971, pp 215-225. 

7n f^i 2 

^ = [log if) + log (-L) + log (-^) + 1 .5]^ + (7) 

* U* 

-4 -1 

where f = coriolis parameter ~ 10 s 
6 - geostrophic wind (m/s) 

Equations 4, 6 and 7 can be used to obtain radar reflectivity as a function 
of the surface layer sensible (H) and latent heat flux (E). Preliminary 
results are shown in Figure 1 for 50 and 300 m above ground. Reflectivity is 
plotted as a function of Bowen ratio (H/E) for curves of constant available 
energy (H+E). The values of available energy are typical for a clear spring 
day and yield reasonable magnitudes of reflectivity from clear air returns. 

As Bowen ratio increases (increase in H, decrease in E associated with soil dry- 
ing) there is a large drop in reflectivity. This trend continues until H be- 
comes approximately twice E. Reflectivity then increases if Bowen ratio becomes 
any higher. The minimum in reflectivity is a result of the negative effect of 
temperature-water vapc • variance in equation 2. 

Conclusions: The reflectivity (n) appears quite sensitive to: 

1) change in Bowen ratio (H/E) when available energy (H+E) is constant. 

2) change in available energy when Bowen ratio is constant. 
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The first case may occur when cloud cover and surface albedo are nearly 
uniform over an area with large differences in soil moisture content. (Soil 
moisture differences may be due to previous rainfall distribution, vegetation, 
and soil conditions.) The available energy can be estimated from the net 
radiational flux (see equation 8). 


R + F = H + E (8a) 

where R = Net radiational flux at earth's surface 

F = Heat flux through conduction from subsurface (F«R during 
periods of high solar radiation) 

R = Solar (shortwave) flux • (1 - albedo) (8b) 

+ longwave flux from atmosphere 
- longwave flux emitted from earth's surface 

Given R, radar reflectivity could be used to map Bowen ratio over the area. 

In other words, both evaporation and heating rate could be obtained with the 
use of radar! 

The second case could occur over an area of nearly uniform soil moisture 
but with differences in albedo or cloud cover. For this condition, both the 
heating rate and evaporation could again be determined from radar reflectivity 
if the Bowen ratio is measured or estimated. 
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STATISTICAL CONSIDERATIONS IN THE ESTIMATION OF WIND FIELDS 
FROM SINGLE DOPPLER RADAR AND APPLICATION 
TO PRE-STORM BOUNDARY LAYER OBSERVATIONS 

Albert 0. Koscielny 

SUMMARY 

Various methods have been suggested for the retrieval of the vector wind 
field from radial velocity data. In all of these methods assumptions must be made 
about either the temporal or spatial structure of the vector wind field. The most 
common such assumption is that the vector wind field is spatially linear and time 
invariant. This assumption has been used for the VAD (Velocity Azimuth Display) 
method of Browning and Wexler (1968), the VARD (Velocity ARea Display) of Easterbrook 
(1975), and the VVP (Velocity Volume Processing) of Waldteufel and Corbin (1979). 

An attractive feature of these techniques is that direct estimates of the 
kinematic properties of the mesoscale wind field (i.e., divergence, deformation, 
vertical shear) can be made from a single Doppler radar. Because these properties 
are important for thunderstorm development, their accurate estimation would be 
a powerful tool for the analysis of the pre-storm environment. 

Waldteufel and Corbin (1979) investigate the application of the VVP to large 
volumes, i.e., a full 360® of azimuthal scan. They find that on this scale, the 
assumption of linearity is a basic limitation to the application of the VVP. In 
this report, only sectors on the order of 30° azimuthal extent are used. It is 
assumed that linearity is a reasonably good approximation for scales on the order 
of 30 km, especially in pre-storm environments. The analysis of radial velocities 
is not straightforward because of the dependence of the accuracy of the estimates 
on the model proposed for the radial velocities and on the geometry of the analysis 
volume. Several methods have been examined and were found to produce estimates that 


1-10 



were biased (VARD) or whose variances were too large (VVP). 

The well known statistical regression theory was used to show that the analysis 
of single Doppler velocities from small volumes is not straightforward but must be 
tailored to specific applications. Considerations from regression theory were used 
to design a model and an analysis volume, termed a modified VVP, that allows the 
estimation of low level divergence with an accuracy of about 3xl0~^s"^ from actual 
radar data. 

The modified VVP was applied to a pre-storm data set for 1530-1630 CST on 
June 19, 1980. The divergence fields from this analysis were found to be reasonable 
mesoscale patterns (See Fig. 1). The fields derived independently from the Norman 
and Cimarron Doppler radars agreed fairly well. Lastly, the areas of convergence 
are areas where cumulus clouds or thunderstorms later develop. 
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Figure 1. Divergence field produced by VVP of single Doppler data 

acquired in a pre-storm environment. Elevation angle ==0.6°; range rings 
are 20 km apart. Divergence is multiplied by 10^. 
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APPENDIX J 


REPORT BY OR. L. RANDALL KOENIG, 
NSF, Washington, D. C., 
(Convective Phenomena Team Member) 


NATIONAL SCIENCE FOUNDATION 
WASHINGTON D C 20550 

DIVISION OF ATMOSPHERIC SCIENCES 


September 2, 1980 


Dr. George H. Fichtl 
ES85 

Marshall Space Flight Center 
Huntsville, Alabama 35812 


Dear George: 


Enclosed are some notes on three observations that could be done using 
the Doppler Lidar. One, the cooperative effort at the BAO could be 
used in proving the system - finding its accuracy and limitations. 

The second, the description of the wind flow around St. Louis probably 
is not appropriate for the FY 81 studies but I would give it consider- 
ation for later studies. If your system operates well and a reasonably 
complete data set on wind flow in St. Louis in spring or summer is 
obtained, many people will be interested. The last suggestion, the 
cirrus cloud study might be useful as a target of opportunity type of 
study. The results could be interesting or they may be uninteresting; 
perhaps a look at cirrus data gathered in your CAT program would be 
useful in developing this study. 

I agree that the anvil mass buduet studies and the lower level flows 
around large convective storms will make interesting studies and I 
suppose someone will write them up. 


I hope these notes are of use to you. 

Sincerely yours. 


L. Randall Koenig I 
Associate Program Director 
for Meteorology 
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1. COOPERATIVE 


Purpose: 


Rationale: 


Procedure: 


EFFORT AT BAO 


To check out instrumentation and contribute to 
boundary layer investigations at BAO. 


Near Boulder a high, well -instrumented tower 
exists, the Boulder Atmospheric Observatory (BAO). 
The tower, with instrumentation at various levels 
is used for boundary layer observations and 
scientists attempt to use these data to develop 
boundary layer models, interpret three-dimensional 
wind fields, etc. This instrumentation should be 
useful for verifying the accuracy of the Doppler 
Lidar, and if sufficiently accurate, the Lidar 
data should be very helpful to the scientist using 
the BAO for establishing the three-dimensional flow. 


Coordination with users of the BAO is required. 

(BAO is operated by the Wave Propagation Laboratory 
of NOAA at Boulder; contact is Dr. J.C. Kaimal 
FTS 320-6261.) A step-wise altitude ascent to 
gather the horizontal winds as a function of alti- 
tude is required (altitudes of data gathering should 
be those at which instruments are placed). 
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2. UNINTENDED WEATHER MODIFICATION 


Purpose: 


Hypothesis: 


What to look 
for: 


Where : 


Procedure: 


To detect mechanisms by which regions of industry 
and urbanization modify weather. 


Industry and urbanization modify weather by 
changing the surface roughness, albedo, etc., and 
by releasing large quantities of heat into the 
atmosphere. Both broad area-wide and essentially 
point sources exist. Urbanization will influence 
mesoscale wind patterns, causing convergence, etc. 
Industrial heat rejection will cause local convec- 
tive anomalies. Both effects should result in 
unintended weather modification manifested by 
increased summer convection and associated phenomena. 


The Doppler lidar would be used to map the wind field 
over the man-altered terrain. Patterns of wind 
convergence on meso to local scales would be sought. 
Influences of topography, heat sources, surface 
roughness changes would be sought. Aerosol distri- 
butions would also be of interest. 


Ideally one would examine different areas, isolated 
industry that emitts large quantities of heat, 
relatively undustry free and flat urbanized areas, 
etc; however, because of the data on anomalous 
weather patterns that exist as a result of the 
Metromex program, I believe St. Louis would be a 
logical site. This region combines urbanization, 
industry having essentially point sources of high 
waste heat rejection and some interesting topography. 
A three-dimensional wind field from the surface to 
the convective condensation height would be very 
useful in sorting out causes of the anomalies 
that are observed. Such information would also be 
useful to numerical modeling attempts that investi- 
gate this matter. 


Data should be gathered in the spring or summer. 

A box like flight path should be used, perhaps 
several rectangles, to cover the needed area. Data 
should be taken at aoproximately 250 m intervals 
from as close to the ground as possible to cloud 
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base height 1000 m (obtain data at four or more 
levels with altitude separation between levels 
increasing with height. Data throughout a day 
would be quite interesting but I suspect that would 
be impossible at the time of the data gathering, 
clouds are not necessary but a sounding would be 
useful . 

The aircraft should carry an IR sensor to map the 
distribution of surface temperature. Other measure- 
ments would also be useful and if you decide to 
follow this suggestion, you might contact Dick 
Dirks at NSF. He may be interested in stimulating 
some additional research. 
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3. CIRRUS CLOUD STUDY 


Purpose: 


What to look 
for: 


When to do the 
experiment: 


Flight procedure: 


Data required: 


To characterize cirrus clouds by their Lidar 
signal. To compare Lidar and visually observed 
characteristics. 


a. Particle distribution in space how homo- 
geneous are the clouds. 

b. Motion field around cirrus (and without it). 

c. Phase change (liquid to ice) in cirrus. Can 
you observe phase transitions? If so, do 
they occur more or less continuously implying 
nuclei are continually being brought into 
the cloud; or that there is a time constant 
(variable) for different nuclei; or that 
some other than heterogeneous nucleation 
process is taking place. 

d. Turbulence within and around cirrus. 

e. Vertical distribution of particles. In 
particular, how do Lidar and visual obser- 
vations differ in this respect? Is one more 
sensitive than the other? 

f. In cirrus that is precipitating ice that 
then evaporates, what is the humidity and 
temperature in the region of evaporation? 

This information would be useful for modeling 
the cirrus generating (seeding) cell concept. 


Target of opportunity. Check out phase around 
San Francisco. Periods of going from one site to 
another site. 


Fly at cirrus level - make measurement from the 
summit to the "base" in steps that seem appro- 
priate depending upon the thickness of the cloud 
(steps of several hundred meters I suspect would 
be appropriate). 


Standard. If it becomes clear that the motion 
field in cirrus is undetectable by the system, 
one might elect to scan only in one direction to 
measure solely the intensity and therefore parti- 
cle distribution in space. 
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